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ABSTRACT 

We present the first interferometric polarization maps of the NGC 2024 FIR 5 molecular core obtained 
with the BIMA array at approximately 2" resolution. We measure an average position angle of — 60°±6° 
in the main core of FIR 5 and 54°±9° in the eastern wing of FIR 5. The morphology of the polarization 
angles in the main core of FIR 5 suggests that the field lines are parabolic with a symmetry axis approx- 
imately parallel to the major axis of the putative disk in FIR 5, which is consistent with the theoretical 
scenario that the gravitational collapse pulled the field lines into an hour-glass shape. The polarization 
percentage decreases toward regions with high intensity and close to the center of the core, suggesting 
that the dust alignment efficiency may decrease at high density. The plane-of-sky field strength can be 
estimated with the modified Chandrasekhar- Fermi formula, and the small dispersion of the polarization 
angles in FIR 5 suggests that the magnetic field is strong (> 2 mG) and perhaps dominates the turbulent 
motions in the core. 

Subject headings: ISM: magnetic fields - ISM: individual: NGC 2024 - polarization- star: formation - 
techniques: interferometric 



1. INTRODUCTION 

Magnetic fields are thought to play a significant role 
in all stages of star formation (e.g., recent reviews by 
Mouschovias & Ciolek 1999; Shu et al. 1999). However, the 
magnetic field is the most poorly measured quantity in the 
star formation process. Observations of the linear polar- 
ization from the thermal emission of magnetically aligned 
dust grains provide a relatively easy approach to explore 
the magnetic field morphology (Heiles et al. 1993). Such 
observations give the field direction in the plane of the 
sky perpendicular to the direction of polarization (Davis 
& Greenstein 1951; Roberge 1996). 

Information on the magnetic field morphology is use- 
ful for testing the predictions of theoretical models and 
simulations. The standard star formation theory predicts 
certain morphological evolution of magnetic fields. The- 
ory predicts that molecular clouds will tend to be flat with 
their minor axes parallel to the field lines, because mag- 
netic fields prevent collapse perpendicular to the direction 
of the field lines (Mouschovias 1976). As contracting cores 
form, the field morphology achieves an "hourglass" shape 
with a collapsing accretion disk ^100 AU at the "pinch" 
and a magnetically supported envelope ~1000 AU (Fiedler 
& Mouschovias 1993; Galli & Shu 1993). Furthermore, the 
rotation of disks may twist field lines into the direction 
along the disk and form a toroidal morphology (Holland 
et al. 1996). Numerical simulations performed by Ostriker, 
Gammie, & Stone (1999) show that the field morphology is 
more random for larger ratios of the thermal to magnetic 
energy. 

Because of the high sensitivity requirement in polariza- 
tion observations, most previous observations of even the 



nearest star-forming regions have been made with single- 
dish telescopes whose large beams cover a region greater 
than the physical extent of protostellar cores. For exam- 
ple, the angular resolution of the James Clerk Maxwell 
Telescope (JCMT) at 850 fim is ~14", which corresponds 
to ^6000 AU at the distance of the Orion Molecular 
Cloud. Therefore, in order to test theoretical models 
and simulations of star formation, it is essential to ob- 
tain high-resolution observations of magnetic fields. The 
best available approach to acquire high resolution infor- 
mation is to conduct interferometric observations. Pio- 
neer interferometric polarization observations at millime- 
ter wavelengths have been done with the Owens Valley 
Radio Observatory (Akeson & Carlstrom 1997). Recently, 
the Berkeley-Illinois-Maryland- Association millimeter ar- 
ray (BIMA) has been able to successfully provide extended 
polarization maps with high resolution up to 2" (Rao et 
al. 1998; Girart, Crutcher, & Rao 1999; Lai et al. 2001, 
hereafter Paper I). 

NGC 2024, a massive star formation region in the Orion 
B Giant Molecular Cloud (distance ~ 415 pc: Anthony- 
Twarog 1982), is a good candidate to explore the magnetic 
field structure at small scales. It contains a luminous H n 
region with a north-south molecular ridge at its center cor- 
responding to the dust lane in the optical image. Mezger 
et al. (1988, 1992) identify seven dust cores in the molec- 
ular ridge at 1300 and 350 /xm and assign designations 
using the acronym 'FIR'. They interpret these cores as 
isothermal protostars in the stage of free-fall contraction. 
However, detailed studies of FIR4/5/6 show that these 
cores contain either a near-infrared source or strong out- 
flows that are traditionally related to protostellar cores 
in a more evolved stage (Moore & Chandler 1989; Chan- 
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dler & Carlstrom 1996). Our interest here focuses on the 
brightest core, FIR 5. With a highly collimated unipo- 
lar molecular outflow extended over ~5' south of the core 
(Richer et al. 1992), FIR 5 appears to be the most evolved 
object among the FIR cores in NGC 2024. Continuum ob- 
servations at 3 mm by Wiesemeyer et al. (1997) resolve 
FIR 5 into two compact cores, and they suggest FIR 5 is a 
binary disk with an envelope. The magnetic field structure 
in NGC 2024 has been mapped with VLA OH and Hi Zee- 
man observations by Crutcher et al. (1999) and with far- 
infrared dust polarization by Dotson at al. (2000). Both 
observations have beam sizes larger than the FIR 5 core. 
Our observations provide new information on detailed field 
morphology in the FIR 5 core. 

2. OBSERVATIONS AND DATA REDUCTION 

The observations were carried out from 1999 March 
to 2001 February using nine BIMA antennas with 1-mm 
Supcrconductor-Insulator-Superconductor (SIS) receivers 
and quarter- wave plates. The digital correlator was set up 
to observe both continuum and the CO J =2^1 line simul- 
taneously The continuum was observed with a 750 MHz 
window centered at 226.9 GHz in the lower sideband and 
a 700 MHz window centered at 230.9 GHz in the upper 
sideband. Strong CO J=2-l line emission was isolated in 
an additional 50 MHz window in the upper sideband. The 
primary beam was ^50" at 1.3 mm wavelength. Data 
were obtained in the D, C, and B array configurations, 
and the projected baseline ranges were 4.5-20, 5-68, and 
5-180 kilowavelengths. The integration time in the D, 
C, and B array was 3.4, 19.1, and 11 hours, respectively. 
The D-array observations were made first on 1999 March 
with a pointing center between NGC 2024 FIR 5 and FIR 6 
(a 20 oo = 5 h 41 m 45!0, 6 2 ooo = -l o 55'55'.'0). Because the 
D-array data showed significant amount of polarized flux 
toward FIR 5 near the edge of the primary beam, the 
follow-up C- and B-array observations were therefore made 
with a pointing center at FIR 5 (a 20 oo = 5 /l 41 m 44!35, 
<5 2 ooo = -l°55'42'/0). 

The BIMA polarimctcr and the calibration procedure 
are described in detail in Paper I (also see Rao et al. 2001, 
in preparation). The average instrumental polarization of 
each antenna was 5.6%±0.4% for our observations. The 
Stokes I image of the continuum was made by mosaicing 
the B-, C-, and D-array data with Briggs' robust weight- 
ing of 0.5 (Briggs 1995; Sault & Killeen 1998) to obtain a 
smaller synthesized beam (HPBW 2'/3xT.'4, PA=9°) with- 
out losing significant amount of flux. D-array data were 
needed to recover the extended emission in order to obtain 
better measurements of the polarization percentage. The 
Stokes Q and U images were made with B- and C-array 
data only, because the precise pattern of the off-axis po- 
larization across the primary beam is unknown. Rao et 
al. (1999) measured the off-axis instrumental polarization 
at selected positions and showed that it only provides an 
uncertainty less than 0.5% (comparable to the uncertainty 
of the leakages) within 4/5 of the primary beam; therefore, 
the off-axis polarization calibration was ignored. Natural 
weighting was used to produce the maps of the Stokes Q 
and U in order to obtain the highest S/N ratio, and the 
resulting synthesized beam was 2'.' '4x1'.' 4 with PA=8°. 

Maps of Stokes Q and U were deconvolved and binned 



to approximately half-beamwidth per pixel (l'.'2x0'.'6) to 
reduce oversampling in our statistics. These maps were 
combined to obtain the linearly polarized intensity (I p ), 
the position angle (0), and the polarization percentage (p), 
along with their uncertainties as described in Section 2 of 
Paper I. When weighted with I p , the average measurement 
uncertainty in the position angle for our observations was 
6?0±2?3. 

3. RESULTS AND ANALYSIS 

Figure 1 shows the B+C-array Stokes I map (con- 
tours) superposed on the D-array Stokes I map (grey- 
scale). The B+C-array map resolves FIR 5 and FIR 6 
into several clumps with the highest resolution ever ob- 
tained (l'/6xT.'2, PA=11°). The clumps with peak flux 
stronger than 7 07 (laj = 3.3 mJybeam -1 ) are named 
and tabulated in Table 1. Figure 2 displays the mosaiced 
map of NGC 2024 FIR 5 overlaid with polarization vectors. 
This map contains B, C, and D array data. Polarization 
vectors are plotted at positions where the observed lin- 
early polarized intensity is greater than 3a7 p (lerj =2.1 
mJybeam -1 ) and the total intensity is greater than 3a j 
(lcr/ = 7.4 mJybeam^ 1 ; note that ctj is dominated by in- 
complete deconvolution rather than thermal noise) . Under 
these criteria, the polarized emission extends over an area 
of ^8 beam sizes. Table 2 lists the polarization measure- 
ments in NGC 2024 FIR 5 at selected positions separated 
by approximately the synthesized beamsize. The distribu- 
tions and the model fitting of the polarization angle are 
plotted in Figure 3 and Figure 4 respectively, and the dis- 
tribution of the polarization percentage is shown in Figure 
5. 

3.1. Continuum and Polarized Emission 

The continuum emission of FIR 5 can be separated 
into three main components - FIR5:Main, FIR5:NE, and 
FIR5:SW (Figure 1). FIR5:Main are resolved into seven 
clumps, and the seven clumps are assigned designations 
with 'LCGR' (Table 1). Only the two brightest clumps 
(FIR5:LCGR 4 and 6) were previously identified at 96 
GHz (Wiesemeyer et al. 1997; note that we have renamed 
these two clumps). Wiesemeyer et al. (1997) suggested 
that FIR 5 is a binary disk; however, our high resolution 
data present complex morphology in FIR 5. Our data also 
show that the continuum of FIR 6 consists of two compact 
sources. The three main components of FIR 5 and the two 
compact sources of FIR 6 are clearly seen in the 450 /im 
continuum maps of Visser et al. (1998), suggesting that 
they are dust dominated sources. 

Figure 2 shows that our detection of the polarized flux in 
FIR 5 is mainly distributed in a ^10"x4" strip along the 
major axis of FIR5:LCGR 4 extending into FIR5:LCGR 
2, which is roughly perpendicular to the disk previously 
proposed by Wiesemeyer et al. (1997). The peak of the 
polarized flux occurs at a position about 1'.'8 north of the 
continuum peak between FIR5:LCGR 3 and 4. Detec- 
tions with I p > 3cr j p were also made toward positions near 
FIR5:LCGR 7 and in FIR5:NE. 

3.2. Polarization Angle Distribution 

The histogram of polarization angles in FIR 5 is shown 
in Figure 3. Most of the polarization detected is associ- 
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ated with FIR5:LCGR 4. Toward this source the average 
position angle is — 60°±6°, and the angles appear to de- 
crease from south to north from —7° to —71°. Compact 
polarization is also observed associated with FIR5:NE and 
FIR5:LCGR 7. Both these sources show a similar po- 
sition angle (average = 54° ±9°) which is very distinct 
from that of the vectors in FIR5:LCGR 4. Our results 
are consistent with recent JCMT observations at 850 /im, 
which also show <fi ~ 50° in FIR5:NE, low polarization 

around FIR5:LCGR 6, and <f> 60° in FIR5:LCGR 4 

(Matthews et al. 2001, private communication). 

Since the magnetic field direction inferred from the 
dust polarization is perpendicular to the polarization 
vectors, the variation of the polarization angles around 
FIR5:LCGR 4 suggests that the magnetic field lines are 
curved. The curved field lines can be successfully modeled 
with a set of parabolas with the same focal point, and 
the best model with minimum \ 2 is presented in Figure 
4. The symmetry axis of the best model is at —77°, which 
is consistent with the position angle of the line connect- 
ing FIR 5: d and f at —67°. The histogram in Figure 4 
shows the distribution of the deviation between our mea- 
surements and the best model, which shows a Gaussian- 
like distribution with a dispersion of <5</> i, s =9?2±3?3. Af- 
ter deconvolving the measurement uncertainty in this re- 
gion (a-0 = 5?8 ± 2?2) from the observed dispersion, we 
obtain the intrinsic dispersion of the polarization angles 
(50 = 7?1±2?5. 

3.3. Polarization Percentage Distribution 

The average polarization percentage is 6.0%±1.5% in 
the main core and 10.5%±4.7% toward the eastern side 
of FIR 5. The polarization percentage distribution is plot- 
ted in Figure 5. These plots show that the polarization 
percentage decreases toward regions of high intensity as 
well as toward the cloud center. Based on the apparent 
close correlation in both Figure 5(a) and 5(b), we per- 
form a least-squares fit on log 10 p versus log 10 / and p 
versus the distance from the peak of FIR5:LCGR 4 (R). 
We obtain the following results: (1) for all data points, 
iogioP = (-1-88 ± 0.01) - (0.83 ± 0.01) x log 10 / with a 
correlation coefficient of —0.86, and (2) in the main core, 
p = (3.60±0.63)%+(1.45±0.16)i? with a correlation coef- 
ficient of 0.97. Because higher intensity and smaller radius 
both imply higher density, our results suggest that the po- 
larization percentage decreases toward the high density re- 
gion. This conclusion is consistent with what we reported 
in Paper I for the W51 el/e2 cores, and our interpretation 
was the decrease of the dust alignment efficiency toward 
high density regions. 

4. DISCUSSION 

4.1. Magnetic Field and Molecular Core Morphology 

In §3.2, we showed that the field lines in NGC 2024 FIR 5 
can be modeled with a set of parabolas which may repre- 
sent part of an hour-glass shape. The reason that the rest 
of the hour-glass shape is not detected could be due to the 
relatively lower dust column density to the east of the re- 
gion where polarization has been detected. The hour-glass 
geometry has been predicted by theoretical work and sim- 
ulations (Fiedler & Mouschovias 1993; Galli & Shu 1993). 
When enough matter collapses along field lines and forms 



a disk-like structure, the gravity in the direction along the 
major axis of the disk would eventually overcome the sup- 
porting forces and pull the field lines into an hour-glass 
shape. Observations of magnetic fields have found sev- 
eral cases consistent with this geometry, such as in W3 
(Roberts et al. 1993; Greaves, Murray, & Holland 1994), 
OMC-1 (Schleuning 1998), and NGC 1333 IRAS4A (Gi- 
rart, Crutcher, & Rao 1999). Our results in NGC 2024 
FIR 5 provide a possible supporting example to the theo- 
retical models at a small scale of few thousand AU. This 
scale is much smaller than the hour-glass shapes observed 
in W3 and OMC-1 (10 4 - 10 5 AU), and is comparable to 
that in NGC 1333 IRAS 4A. 

The symmetry axis of our hour-glass model is approxi- 
mately parallel to the binary disk previously proposed by 
Wiesemeyer et al. (1997). Such a coincidence could be 
an indirect evidence for the existence of a collapsing disk. 
Although our high resolution map shows that FIR5:Main 
is more complicated than a binary disk, it is still possible 
that there is an east-west disk consisting of four clumps 
(FIR5:LCGR 1, 4, 6, and 7), with the remaining three 
clumps either collapsing into the disk or being ejected 
through outflows. FIR5:LCGR 2 is at - -145° direction 
with respect to FIR5:LCGR 4, which is close to the unipo- 
lar outflow at ~ —170°. Since outflows could be deflected 
by magnetic fields away from the protostellar cores (Girart, 
Crutcher, & Rao 1999), it is possible that FIR5:LCGR 2 
traces the direction of the outflow near the central core. 
Further high resolution kinematic study is needed to ex- 
amine the above speculation. 

4.2. Estimation of the Magnetic Field Strength 

Although the magnetic field strength cannot be di- 
rectly inferred from polarization of dust emission, Ostrikcr, 
Stone, & Gammie (2001) show that the Chandrasekhar- 
Fcrmi formula (Chandrasekhar & Fermi 1953) modified 
with a factor of 0.5 can provide accurate estimates of 
the plane-of-sky field strength under strong field cases 
(S(j) < 25°) for their simulations of magnetic turbulent 
clouds. Therefore, the projected field strength (B p ) can 
be expressed as 



B p = 0.5 y^p 
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(1) 

where p is the average density, 5vi os is the rms line-of- 
sight velocity, the number density of molecular hydrogen 
n H2 — 2.33 m H n H , and the linewidth Av = \/81n2 Svi os . 

To estimate B p in NGC 2024 FIR 5, we must first care- 
fully determine the density of the dust core and the tur- 
bulent linewidth. Mezger et al. (1992) derived high den- 
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in FIR 1-7 from dust emission, 



but their value is significantly higher than that derived 
from molecular studies (^10 6 cm~ 3 from CS: Schulz et al. 
1991; ~10 7 cm~ 3 from C 18 0: Wilson, Mehringer, & Dickel 
1995). It is possible that depletion of molecules in the 
dense core causes molecules to only trace the envelope of 
the dense core; on the other hand, Chandler & Carlstrom 
(1996) showed that the kinetic temperature in the FIR 
cores is higher than what Mezger et al. (1992) assume, 
implying a lower density. Mangum, Wootten, & Barsony 
(1999) studied the kinetic temperature in NGC 2024 with 
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multi-line observations of formaldehyde (H2CO), and their 
results supported the arguments of Chandler & Carlstrom 
(1996). We therefore adopt n H2 ~ 2 x 10 6 cm" 3 for FIR 5 
from Mangum, Wootten, & Barsony's large velocity gradi- 
ent (LVG) calculation. We also adopt a weighted lincwidth 
of Aw = 2.04 ± 0.03 km s _1 at the position at the peak of 
FIR 5 from Mangum, Wootten, & Barsony (1999). Note 
that, just as for other molecules, the peak of H 2 CO emis- 
sion does not coincide with the continuum peak; there- 
fore, the H 2 CO emission may come from the envelope of 
the FIR 5. This may be an advantage because the H2CO 
lincwidth is less likely contaminated by the possible dy- 
namical motions in the core and better represents the tur- 
bulent motion. 

In §3.2, we calculated the angle dispersion 5<fi = 7?1±2?5 
after taking out the systematic field structure in order 
to identify the dispersion purely from the Alfvenic mo- 
tion. Along with the parameters discussed in the previous 
paragraph, we obtain B p ~3.5 mG. We can also calculate 
the lower limit of B p using the largest angle dispersion, 
which is the dispersion before taking out the parabolic 
model (S(j) =13?1±4?6); thus, B p>min -1.9 mG. Our esti- 
mate of the planc-of-sky field strength is much larger than 
the line-of-sight field strength measured by Crutcher et al. 
(1999) from OH Zeeman observation, which is ~ 65^iG at 
FIR 5. Except in the unlikely case that the magnetic field 
direction lies almost on the plane of the sky, the low field 
strength detected with the Zeeman measurement can be 
explained by (1) the beam-averaging over small scale field 
structure, and/or (2) that OH line does not trace high 
density region of the dust cores. 

4.3. Turbulent Energy vs. Magnetic Energy 

The small dispersion of the polarization angles observed 
by us may imply that the turbulent motion is not strong 
enough to disturb the magnetic field structure. In order 
to quantitatively discuss the the relative importance of the 
turbulent and magnetic energy in NGC 2024 FIR 5, we cal- 
culate the ratio of the turbulent to magnetic energy 

Pturb = -^r, (2) 

where a tur b is the turbulent linewidth and Va = |B|/^/47rp 
is the Alfven speed. Statistically, < B 2 >= | |B| 2 , thus 
Va can be estimated from Eq. (1). Therefore, 

Pturb (3) 

In the general conditions prevalent in molecular clouds, the 
thermal linewidth a t hermai is much smaller than the tur- 
bulent linewidth, so 5vf os = a 2 turb + a 2 thermal w a 2 turb . The 
ratio of the turbulent to magnetic energy simply depends 
on the polarization angle dispersion, 

/^, & «^ 2 =8xl0- 4 (^) 2 . (4) 

The 0.5 factor in Eq (1) was obtained under the as- 
sumption that the dust alignment efficiency is uniform 
throughout the cloud. If dust alignment efficiency is lower 
in the regions with higher density as we suggest in §3.3, 
Pturb would be larger for the same 8(j> and the correction 
depends on the degree of depolarization. For our case, 
the most conservative estimate of the angle dispersion is 



6<fi < 13? 1, which leads to a small turbulent to magnetic 
energy ratio, fiturb < 0.14. Therefore, the magnetic field 
most likely dominates the turbulent motion in the core 
region of NGC 2024 FIR 5. This is also the case in W51 
el/e2 (see Paper I). 

5. CONCLUSIONS 

High-resolution polarization observations are needed 
to explore the detailed magnetic field structure in star- 
forming cores. We have obtained a continuum map of 
NGC 2024 FIR 5 and FIR 6 and a polarization map of 
NGC 2024 FIR 5 with the highest resolution ever obtained 
(~2"). Our observations resolve FIR 5 and FIR 6 into sev- 
eral continuum clumps. The information revealed by our 
polarization observations of FIR 5 are summarized below: 

• Extended polarization is detected associated with 
the main core at — 60°±6°. Compact polarization 
is also observed toward the eastern side of FIR 5 at 
54°±9°. The polarization is low between these two 
regions and the secondary intensity peak of FIR 5 
has an upper limit of 6%. 

• The magnetic field lines in the core are systemati- 
cally curved with a symmetry axis close to the ma- 
jor axis of a putative disk. This is consistent with 
an hour-glass morphology for the magnetic fields 
predicted by theoretical works. 

• The polarization percentage decreases toward re- 
gions with high intensity and short distance to the 
center of the core. The tight correlations imply 
that the depolarization is a global effect and may 
be caused by the decrease of the dust alignment 
efficiency in high density regions. 

• The small dispersion of the polarization angles in 
the core suggests that the magnetic field is strong 
(> 2 mG) and the ratio of the turbulent to magnetic 
energy is small (flturb < 0.14). Therefore, the mag- 
netic field most likely dominates turbulent motions 
in NGC 2024 FIR 5. 
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Table 1 

Compact continuum sources in FIR 5 and FIR 6 





Ck2000 


^2000 


J a 


P 


Source 


h m s 


o / // 


(mJybeam -1 ) 


(%) 


FIR5:LCGR 1 


5 41 44.00 


-1 55 40.7 


46 


< 17 


FIR5:LCGR 2 


5 41 44.10 


-1 55 44.0 


64 


10.3 


FIR5:LCGR 3 


5 41 44.23 


-1 55 36.8 


39 


< 25 


FIR5:LCGR 4 b 


5 41 44.25 


-1 55 40.8 


293 


3.8 


FIR5:LCGR 5 


5 41 44.32 


-1 55 35.7 


40 


< 21 


FIR5:LCGR 6 c 


5 41 44.48 


-1 55 42.2 


107 


< 6 


FIR5:LCGR 7 


5 41 44.69 


-1 55 43.5 


55 


< 14 


FIR6 c 


5 41 45.13 


-1 56 04.2 


67 




FIR6 n 


5 41 45.17 


-1 56 00.3 


26 





"The peak intensity is measured from the B+C+D-array map. 
6 FIR5-w in Wiesemeyer et al. (1997) 
c FIR5-e in Wiesemeyer et al. (1997) 



Table 2 

Polarization Measurements in NGC 2024 FIR 5 



Position a 


Stokes / b 


Polarization 


Polarization 


Note 


(",") 


(mJybeam -1 ) 


Percentage(%) 


Angle (°) 




(-0.3, 4.5) 


26.6±7.3 


23.1±9.5 


-72±9 




(-1.1, 3.0) 


132.8±7.3 


18.7±1.9 


-71±2 


Polarization Peak 


(-1.2, 1.2) 


300.6±7.3 


3.9±0.7 


-64±5 


Intensity Peak 


(-1.8, 4.5) 


41.9±7.3 


19.5±6.1 


-66±7 




(-2.6, 1.2) 


153.7±7.3 


4.9±1.4 


-39±8 




(-2.6,-0.8) 


137.1±7.3 


6.5±1.6 


-39±7 




(-2.6,-2.7) 


45.6±7.3 


15.4±5.2 


- 8±9 




(-3.9,-0.9) 


53.2±7.3 


20.0±4.8 


-42±6 




(-4.8,-2.4) 


39.1±7.3 


17.3±6.3 


-39±9 




(10.2, 8.9) 


31.7±7.3 


27.1±9.2 


63±7 


FIR5:NE 


( 4.1,-0.5) 


75.7±7.3 


9.3±2.9 


31±9 


near FIR5:LCGR 7 






6.0±1.5 


-60±6 


average of the main core 






10.5±4.7 


54±9 


average of the eastern side of FIR 5 



"Offsets are measured with respect to the phase center: Q 2 000=19 ,l 23 m 44 , ;2, 5 2 ooo=14 30'33'.'4. 
6 Stokes I is measured from the B+C+D-array map. 
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a (J2000) 

Fig. 1. — Stokes I maps of NGC 2024 FIR 5 and FIR 6. The contours show the total intensity (Stokes I) in the B+C-array data at 
—5,-3,3,5,8,12,17,23,30,40,50,60 a levels (lcr = 3.3 mjy beam -1 ). The color of the contours is chosen to show better contrast on the 
background. The field of view is represented by the northern dotted circle. The grey-scale shows the D-array data and the field of view is the 
southern dotted circle. The grey ellipses in the upper left corner are the synthesized beams for these two sets of data, which are l'.'6x 1'.'2 with 
PA=11° for the B+C-array map and ll"x8" with PA=16° for the D-array map. The triangles mark the positions of compact continuum 
sources identified in the B+C-array map (Table 1). 
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05 h 41 m 45 s 5 45 s 44 s 5 44 s 43 s 5 

a (J2000) 

Fig. 2. — Polarization map of NGC 2024 FIR 5. The contours represent the mosaiced Stokes I map at —5, —3, 3, 5, 10, 20, 30, 40<r levels. 
The 1 a noise level of Stokes / is 7.3 mjybcam~ 1 and the beam HPBW is 2'.'3xl'.'4 with PA=9°. The line segments are polarization vectors 
(E- vectors), and their lengths are proportional to the polarization percentage with a scale of 10% per arcsec length. The triangles mark the 
positions of compact continuum sources listed in Table 1. 
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Fig. 3. — The distribution of the position angles of the polarization vectors in NGC 2024 FIR 5. The vertical axis is the total number of 
the measurements in the bin. 
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Fig. 4. — Model fitting of the magnetic field morphology using a set of parabolas with the same focus point. The line segments are 
polarization vectors which are perpendicular to the magnetic field directions. The dotted lines shows the least-squares model, the dashed 
line is the symmetry axis, and the star is the common focus point. The histogram shows the distribution of the deviation between the field 
directions and the model. 
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Fig. 5. — Polarization percentage distribution. The filled dots represent the data associated with the main core of FIR 5, and the open dots 
represent the data associated with the eastern side of FIR 5. The dotted lines are the best least-squares fits for the set of data, (a) shows the 
polarization percentage vs. the total intensity. The error bar of each data point is the measurement uncertainty. We plot the total intensity 
increasing to the left to emphasize that high intensity corresponds to short distance to the center of the core, (b) shows the polarization 
percentage vs. the distance from the FIR5:LCGR 4 peak. Data with R greater than 6" are not shown, because the average polarized intensity 
is lower than 1.5o\r . The error bar in radius shows the range over which the data are averaged and the error bar in p is the standard deviation 
of the data in the range. 



